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Abstract 

Recently, it has been observed that the effective dipolar interac- 
tions between nuclear spins of spin-carrying molecules of a gas in a 
closed nano-cavities are independent of the spacing between all spins. 
We derive exact time-dependent polarization for all spins in spin- ^ 
ensemble with spatially independent effective dipolar interactions. If 
the initial polarization is on a single (first) spin, Pi (0) = 1 then the 
exact spin dynamics of the model is shown to exhibit a periodical 
short pulses of the polarization of the first spin, the effect being typ- 
ical of the systems having a large number, N , of spins. If N S> 1 , 
then within the period An/g (2Tr/g) for odd (even) iV-spin clusters, 
with g standing for spin coupling, the polarization of spin 1 switches 
quickly from unity to the time independent value, 1/3 , over the time 
interval about (g^/N)^ 1 , thus, almost all the time, the spin 1 spends 
in the time independent condition P\(t) = 1/3 . The period and the 
width of the pulses determine the volume and the form-factor of the 
ellipsoidal cavity. The formalism is adopted to the case of time vary- 
ing nan o- fluctuations of the volume of the cavitation nano-bubbles. 
If the volume V(t) is varied by the Gaussian- in-time random noise 
then the envelope of the polarization peaks goes irreversibly to 1/3 . 
The polarization dynamics of the single spin exhibits the Gaussian ( or 
exponential ) time dependence when the correlation time of the fluctu- 
ations of the nano- volume is larger ( or smaller ) than the {(5g) 2 )~ 1 / 2 , 
where the ((5g) 2 } is the variance of the g(V(t)) coupling. Finally, we 
report the exact calculations of the NMR line shape for the ./V -spin 
gaseous aggregate. 

PACS numbers: 05.30.-d, 76.20. +q 
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1 Introduction 



The nature of ergodicity being of fundamental importance for consonant 
description of statistical mechanics is currently being discussed in the NMR 
context pQ. Spin dynamics is ergodic if the initial polarization prepared at 
a single (first) spin is spread over the system leading, as time proceeds, to 
the spatially uniform distribution of the polarization, as it is expected on the 
basis of a simple physical intuition. Instead, non-ergodic behavior that have 
recently been observed numerically in the nuclear spin- | ID chains under 
the general XYZ spin Hamiltonian |2J enters in the way that the time average 
polarization of the first spin turns out to be several times larger as compared 
to that of any other spin in the chain. This observation of nonergodicity has 
been extended to ID chains and rings under XY Hamiltonian |3] showing 
analytically that the time average polarization of the first spin differs by the 
factor 1.5 -f- 2 from the time average polarization of all other spins in the 
chain. These considerations in ID spin clusters address the problem of the 
nature of the ergodicity to different spin Hamiltonians. Motivated in studying 
of non-ergodic spin dynamics and due to the fact that exact solution is the 
lucky exception in statistical mechanics, our assumption in this paper is that 
spin interactions are regarded to be independent of the spacing between the 
spins rather than having r -3 dependence. 

Recently, the spin Hamiltonian with space-independent spin couplings has 
been applied for exploring the NMR spectra of a gas of spin-caring molecules 
undergoing a fast thermal motion within the non-spherical cavities [4. In 
that report, the authors have arrived at the space- independent effective spin 
couplings by motionally average the exact dipolar Hamiltonian over a uni- 
formly distributed spins' spatial coordinates in a nanometer sized cavities. 
This technique is expected to have a promising application for determining 
the pores' shapes and sizes j3] by the NMR spectra. 

With regard to the effective nuclear spin Hamiltonian with infinite range 
couplings it is noteworthy that this type of interactions has also been pro- 
posed in the theory of nano-electrodes jB] , [7] . Infinite range dipolar nuclear 
interactions are induced indirectly due to the fast energy transfer between 
electron and nuclear spins. On the coarse grained time scale of the fast 
electron spin dynamics, the slow nuclear spin dynamics is governed by an 
effective nuclear spin Hamiltonian with an infinite range interaction cou- 
pling. Quite apart from its importance as a physical model in the NMR 
experiments for a many-spin aggregate in a confined volume [I], [E], [7] and 
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few proton molecules jSj, the model with infinite range spin interactions is 
of a fundamental interest in its own right since this model allows to treat 
3-dimensional case exactly without any reference to an ID spin ordering. 
It represents the quantum non-equilibrium version of the exactly solvable 
equilibrium spin model |DJ , has the mapping to the BCS pairing Hamiltonian 
of the superconductivity [TU], and has long been provided to be the test for 
many-body problems in high spatial dimensions, D ^> 1 . The objective 
of the paper is to present the exact solution of non-ergodic dynamics with 
infinite range spin Hamiltonian in N- spin- | ensemble. 

To our knowledge, the only result reported on this model is that of Waugh 
[TT] who announced (without proof) that the time average polarization of 
the first spin equals exactly and the polarization of any other spin 

is exactly ^ for odd numbered, N , spin cluster. To clarify the problem 
of spin dynamics, the present paper reports a detailed analytical theory of 
the average polarization both for odd and even numbered spin clusters as 
well as it gives the theory of spin dynamics that is entirely missing from 
11J. A condensed form of this paper have been published in [12] • Brief 
overview of the paper is as following. In section 2 , we construct the effective 
nuclear spin Hamiltonian of spin-carrying molecules of the perfect gas in 
the nano-cavity. Section 3 gives formalism required to obtain the exact 
time-dependent polarization. This is followed by Section 4 that discusses 
three issues of the polarization dynamics that amenable to the techniques 
of the Section 2 : firstly, the non-ergodicity of the polarization dynamics of 
a single spin in the nano-cavity; secondly, the polarization dynamics of a 
single spin within a fluctuating nano-bubbles; thirdly, spectral line shape of 
the nuclear spin ensemble. Finally, Section 5 summarizes the conclusions 
of the calculations and confronts the results obtained with known analytical 
results for the XY Hamiltonian. 

2 Effective Nuclear Spin Hamiltonian in Nano- 
Cavity 

The purpose of this section is to construct an effective spin Hamiltonian, 
H e ff , that governs spin dynamics of spin-carrying molecules in a nano-sized 
cavity on the coarse grained temporary scale of the order 10 picoseconds. 
On these space-time scales the effective spin Hamiltonian differs from the 
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exact dipolar Hamiltonian, in particular, the H e ff has a high symmetry that 
permits the exact solution for the spectrum of the many-body spin Hamilto- 
nian H e ff and, as a result, the exact derivation of the polarization dynamics 
of the gas within the nano-cavity. In this section we summarize the main 
ideas of the report [3] , however in derivation of the effective spin Hamiltonian 
H e ff by averaging over spins' spatial coordinates, we generalize the effective 
spin coupling to the case of non-perfect gas in the nano-cavity 

The starting point of the derivation of the operator H e ff is the expression 
for the density matrix 

p(t,{I n ,f n (t),p n (t)}n=i) = ?7(t)p(0,{4f B (0),^(0)}2 =1 )» r - 1 (*) (!) 

for completely specified coordinates {r n (t)}^ =1 and momenta {p n (f)}n=i of 
N spin-carrying molecules. The propagator U (t) is associated with the time 
dependent exact dipolar Hamiltonian (in frequency units) 

N 

H{t) = £ hij(t), tnj(t) = 1 2 %P 2 {cose ij {t))r^{t){I i I j - 3I iz I jz ), (2) 
i<i<j 

where, 7 stands for the gyromagnetic ratio, I na (a = x,y,z) specify the 
spin- I operators, the is the instant polar angle between the vector 

Tij(t) from Ti(t) to Tj(t) and the external magnetic field B. 

A cornerstone fact for a construction of an effective spin Hamiltonian is an 
essential difference between the time scale of the relaxation in the phase space 
r N —p N and the time scale of the spin dynamics under the Hamiltonian of Eq. 
((21). Actually, for the hydrogen gas at room temperature and atmospheric 
pressure, the following estimations hold. The average concentration n fa 
2.7 • 10 19 molecules/cm 3 , the mean free path A = ( 10 cm for 

the radius of a molecule a w 1CT 8 cm and the thermal velocity v w 10 5 
cm/s. Then, a simple order-of- magnitude calculations leads us to expect 
that, for the gas in the cavity of the size £ fa 10 nm, the diffusion coefficient 
T> ~ vi ~ 10 _1 cm 2 /s, the characteristic time scale of the spatial relaxation 
of the gas is ~ C 2 /T> fa 10 -11 s and the characteristic time scale of the 
momenta relaxation towards the Maxwellian distribution is i v ~ X/v ~ 10 -11 
s. These time scales t v and td\f are well separated from the NMR time scale 
^nmr = 10~ 4 ~v~ 10~ 3 s associated with the dipolar interaction in Eq. ©. The 
smallness of the parameter 

e = — — = 10 -7 1, where t re \ = max(t v , t dif ) (3) 
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allows for determining the average nuclear spin Hamiltonian governing the 
behavior of the nuclear spins over a coarse-grained time intervals At obeying 



*rel < At < t nmr . (4) 

The averaging of the exact Hamiltonian over time At is performed to the 
zeroth order in the perturbation expansion in powers of the parameter e 
yielding the average (or effective) Hamiltonian [T5] . 

1 /-a* 

hj = ic7 / (5) 



At Jo 

with the corrections being of the order 0(e l ) . 

The decisive point of the following treatment is the replacement of the 
temporary integration in Eq. with the integration over the spatial coor- 
dinates within the confined region. Equating the temporary averaging with 
the spatial averaging makes sense under the ergodic hypotheses 



6t(dr N ,dp N ) = z -i e -E/kT dr N dp N^ (6) 
where the notations imply that a representative point living in the whole 
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phase space r — p , while moving over the time t , t re \ < f < („ 
spends within the volume dr N dp N only a fraction Z~ 1 e~ E / hT dr N dp N of 
the whole time t , the E being the total energy. The relation (JHJ) incorporates 
the Gibbsian stochastic level of description into the dynamical treatment of 
the molecular collisions. The phase space degrees of freedom are assumed to 
relax to their equilibrium distribution at a given temperature T . 

Introduce equilibrium pair distribution function for molecules 1 and 2 , 

n , ] = J v d 3 r 3 ---J v d 3 r N exp(-U(r N )/kT) 
2[n,r2) J v cPrv • • J v d s r N eM-U(r N )/kTY 1 ) 

where U(r N ) stands for inter-molecular electrostatic interactions, (recall 
that ||[/||/||If || w 10 7 [T5]). Then, taking for granted the ergodicity (0), 
the evolution of spin degrees of freedom on the coarse grained time scale At 
of Eq. (@J) is governed by the static (time independent) effective Hamiltonian 

TV N 

H = Kj = 9 ~ siizijz) (8) 

l<i<j l<i<j 
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with spacing independent pair couplings g for any pair of spins i and j , 

9 = l 2 hj v J y d\i d 3 rj D 2 (r t , r j )P 2 (cos %)r,r 3 . (9) 

The effective operator H involves only the (slow) spin operators, whereas 
the (fast) spatial coordinates of the nucleus ( labeled by indices % and j ) are 
integrated out. On the coarse grained scale At , any nuclear spin "feels" the 
field that is independent of the spatial coordinates of all other spins flying 
within the nano-cavity but depend on the quantum states of that spins. 

The effective spin coupling g encodes an information about the shape 
and the size of the nano-cavity. The primary objective of the proceeding 
discussion is to present the expression for the coupling g of Eq. Q for an 
ellipsoidal nano-cavity. 

For a perfect gas within the nanocavity, the pair distribution function 
equals D 2 (ri,rj) = V~ 2 for spins i and j lying in the cavity r^Tj e V , 
hence the averaging in Eq. Q gives 

g = j 2 hF/V, F = JJ y ^f^P 2 (cos6 12 )r^. (10) 

In this report, use is made of the remarkable fact jl] that the volume, V , of 
the nano-cavity enters into the expression for the effective coupling g of Eq. 
(I10J) . which itself enters into the polarization (defined below in the section 
EJ) giving rise to the dependence of the polarization on the volume of the 
nano-cavity by no means as trivial as merely proportional to the volume. 

The transformation of the coordinates r\ and r 2 to the relative coordi- 
nate r 12 = r\—r 2 and the coordinate of the center of gravity, r = (ri+r 2 )/2 , 
reduces the form-factor F in Eq. (jlOj) to the form 

F = [ d 3 r 12 P 2 (cos6 12 )r^. (11) 

It is convenient to assume that the initial point of the vector T\ 2 = n — 
r 2 starts at the origin of the frame of reference xyz connected with the 
ellipsoid, see Figure 1 . In order to calculate the integral in Eq. (jTTj) in 
the ellipsoid frame of reference, we define by z , ri2 and Ziri2 the polar and 
the azimuthal angles of the vector r\ 2 in the ellipsoid frame of reference 
xyz , respectively, and the angles #z,ri 2 > 0z,r 12 wm be referred to as the 
polar and the azimuthal angles for the vector r 12 in the laboratory frame 
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of reference XY Z . By the " addition" theorem of spherical harmonics, we 
can write for the factor P 2 (cos$i 2 ) in Eq. ([lip the expression P 2 (cos0i 2 ) = 
P 2 (cos 9z,r 12 ) — p2(cos6 l Zj z)P2(cos6 l Ziri2 ) + ... where the dots denote the 
terms depending on the azimuthal angles as exp{iw(0 ZiI . 12 — </> z ,z)} with 
m^O, hence, all these terms vanish after integration over the angle (j) z ,r 12 ■ 
Defining shorthand (ft = z , ri2 , 9 = 9 z>ri2 , a = 9 z Z , we get 

F = P 2 (cosa)F', F' = d<p sin 9 d9P 2 ( cos 9) r^dr 12 . (12) 

Jo Jo Jo 

The Eq. ()12|) shows that the dependence of the form-factor on the angle a 
is factored out. The remaining integral, F' , in Eq. ()12j) coincides exactly 
with the form-factor of the ellipsoid having the z axis along the Z axis of 
the laboratory frame of reference. The radius R(9) stands for the surface of 
the ellipsoidal cavity. To find the integral in Eq. (|12jl . we parameterize the 
ellipsoid {z/a) 2 + (x 2 + y 2 )/b 2 = 1 with the spherical coordinates p and 9 
as follows, p = 6(1 — e 2 cos 2 #) -1 / 2 , where e 2 = 1 — (b/a) 2 . In the spherical 
coordinates, the form-factor F' becomes 



F' = d(j> sin0d0P 2 (cos0)ln 
Jo Jo 



b/ Pmin 



1 -e 2 cos 2 9y/ 2 



, (13) 



here, the auxiliary lower limit P m i n — > is introduced for the convergence 
of the integral. The constant P m i n vanishes from the final form-factor F' . 
Indeed, invoking the definition P 2 (cos0) = (3cos 2 — l)/2 and using the 
integral j\P 2 (x)dx = , we get rid of the dependence on P min in the 
integral F' of Eq. ()13p. thus, reducing the integral F' to the form 

F' = -- f 277 d<f> r sin 9d9P 2 (cos 9) ln(l - e 2 cos 2 9). (14) 
2 Jo Jo 

Straightforward integrating over 9 by parts and combining the result with 
Eq. (JT2~j) yields the sought for form-factor 

f | + 2 (\ - I) fl - i arcthe) , for a > b 
P = PrP 2 (cosa), 1=1 I Y\ \\ % '\ " (15) 

1 |-2(^ + l) (l-iarctg|e|), for a < b . 



For a > b , the e = 1 and 1 = 2/3. For a = 6 , the e -> and Eq. (fT5j) 
gives 7 = 0. For a <^ b , the |e| — > 00 and 7 = | — 2 = —4/3. These 
limiting cases confirm the result reported in jl]. 
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3 Polarization Dynamics 



We regard the spin Hamiltonian H of a iV -spin cluster in a uniform external 
magnetic field B parallel to the Z-axis of a fixed frame of reference XYZ 
and the spatially independent spin couplings g of Eq. (fTO]) . 

— N 9 N 

H iO ^ ' I nz ~t~ ~ ^ ] {C Imzlnz Imxlnx I-mylny} i (^) 
n=l m^n 

here, u; = 7-B denotes the Zeeman frequency and the ( is an arbitrary 
factor. 

The standard way of approaching iV-spin aggregate is to find the polar- 
ization on the n-th spin at a time instant t given initial polarization on the 
1-st spin 

P n (t) = tr {e mt I lz e~ m I nz } /tr {l 2 lz } . (17) 

The Hamiltonian of Eq. (J16)) can be cast into the form (up to the constant 
gN(l - §)/4 ) 

H = uI z + 9 -{C + l)I 2 z - g -l\ (18) 

where the total spin / = J2n=i In and its projection onto the z-axis reads 
Iz = J2n=i Inz ■ The polarization P n (t) of Eq. (fTTjl is unchanged if we modify 
the Hamiltonian H of Eq. |TH|) to the effective one 

H> = - 9 -l\ (19) 

In addition, the equivalence of polarizations P n {t) of all spins, except the 
first spin, as well as the conservation with the time of the total polarization 
J2n=i Pn{t) allow to concentrate the description on the first spin only, 

P 1 (r) = tr {e^ 12 I lz e- lTl ' 2 I lz } /tr {l 2 lz } , (20) 

here, the dimensionless time scale is defined according to r = ^gt . A pow- 
erful tool for investigating the problem of Eq. (J2*U|) is the theory of coupling 
of the angular momenta jl(i| . |17j . To this end, think of the total spin cluster 
as composed of two subsystems A and B . The subsystem A has only the 
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spin I\ = I a and the remaining fragment B of the spin cluster has the spin 
I B , with the total spin being I = Ia + I B . The states of the two subsystems 
A and B are coupled together within the state of the whole system A © B 
through the Clebsch-Gordan (C-G) coefficients, 

\I A ,lB,I,m)= C lZnA;lB,m B \ I AmA)\lBm B ), (21) 

m A=±? 

where Ia = \ and uia = ±2 are ^ ne s P^ n an< ^ ^ s m agnetic quantum numbers 
for the 1-st spin, respectively, and I B and tub — m — rriA stand for spin and 
its magnetic quantum numbers of the fragment B . For Ib = only I = \ 
and m = ±| are allowed. For I B >\, the allowed I and m are as follows, 
I = I B ± i , —J < m < I . The C-G coefficients are, see e.g. [TH] . 

= c?;> m+i = [(/ fl + i +m )/ ( 2/ B + i)]^, 

C^fc^J = C^ m _,=[(/ B + i-m)/(2/ B + l)]^ (22) 

The two pairs of independent variables (I B ,m B ) and (I a = \,m A = ±|) 
will be used for determining the trace in Eq. (J2(Jj) for the whole N spin 

system A © B , 



N B /2 Ib + \ I 

tr{...}= £ tu(I B ) £ E (I A jB,I,m\...\I A ,I B J 1 m) 1 (23) 

/ B =/gi" r=|/ B _i|m=-/ 



where Nb = N — 1 is the number of spins in the fragment B , the minimal 
value of Ib is I™" — for even Nb and I B m = \ for odd Nb ■ The factor 

2Ib+1 I N B + 1 \ 

is the number of ways of grouping Nb individual spin- | into the total spin 
I B • The w(I B ) obeys the property [TH], [IZj, [E| 

V B /2 

(25) 



E ! 

>lmnl \2 



The right-hand side of Eq. ()25|) is the number of states for each allowed 
eigenvalue m B of the fragment B. In order to deal with the diagonal evo- 
lution matrices within the Eq. (J2(Jj) . we introduce additional basis of the 



9 



bra, (Ia, Ib, I', m '\ , and the ket, | I A , Ib, I', m ') , vectors of the Hilbert space 
'H(Ia) ® T~C(Ib) f° r fixed values Ia = 1/2 and the Ib , then make use of the 
completeness of the 2(21 b + 1) orthonormal basis vectors belonging to the 
space H(I A ) ® H(I B ) 



i B +k 



i 



l H{i A )®H(i B ) = E \I A ,I B ,I\m'){I A ,I B ,I',m'\, (26) 

r=\i B -l\m'=-i 



and, finally, insert the representation of the unity of Eq. (|26)l in front of the 
right-most operator I\ z of Eq. (|2*nj) having the matrix elements 

(I A ,I B ,I'm'\Ii z \I A ,lB,I,rn) = 5 mm > ■m A C 1 1 ,m , . C[' m T .(27) 

N a > ' lz| 131 ' ' in, in j-i ± jm , A -I Bj rn'—m A ±,m A ;I B ,m—m A 

m A =±\ 

With these algebraic steps, we immediately get the polarization P\(t) in 
terms of the Clebsch-Gordan coefficients as 

N B /2 

Ib=IT \i B -K<*<'b+\ J£%, 

\I B -^\<I'<I B +h ~ ~ 
xe ir{i(/ + i)-/'(i'+i)}/ y m ^{V , C[ m T )\ (28) 

V * — ' 2> mA > B ' m mA 2' mA < B ' m ~ mA '' 

m A =±\ 

For the term = , only the single pair (I — ^,1' — |) is allowed in 
the sum Eq. (J2"8j) . and for I B > % > the f° ur pairs of (/,/') should be 
distinguished in the sum Eq. ()28|) depending on the sign (+) or (— ) in the 
expressions 

(IJ')=[i = Ib±\,I' = Ib±1). (29) 

Armed with the polarization P\(t) of Eq. (|28jl . we are going to decom- 
pose it into the time-independent part Pi and the oscillating part P° sc (r) , 

P 1 (r) = p- 1 + Pr(r). (30) 

The time-independent contribution Pi to the function P\(t) is provided 
by the quantum numbers m, m' belonging to the states / = I' = Ib ± 5 
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if Ig > k , and by the quantum numbers m, m' belonging to the states 

i = r = \ if i B = o, 

N b /2 1b + \ I 

a = 2-<"»-» Y. E E ( E ^(cj;: AifB , ra _ m J 2 ).(3i) 

/ s =/g m /=|-fs-|l m =- 7 m A =±| 

Our desire now is to sum in Eq. (jHljl over the indexes m and J for the fixed 
value of the Ib ■ To this end, we start with the state Ib = that arises 
for even Nb ( see the comments to Eq. ). For Ib = only I = | is 
allowed, and the partial polarization P\{Ib) in Eq. (J3TJ) reads 

i 

T\(I B = 0) = 2-( iVs - 1 ) W (0) E m 2 . (32) 



Next, we consider the contribution to the P\ of Eq. (|3*T|) from the spin 
Ib > §■ m this situation, I = Ig ± | are allowed and invoking the C-G 
coefficients of Eq. (J22J) . the contributions Pi(Ib) to the Pi of Eq. (|3*Tj) can 
be written in a convenient form as 

T l {I B )=2- N »w{I B ) E (33) 

Combining together the 7\(I B = 0) of Eq. §2$ and the 7\(I B ) of Eq. ® 
results in 

N B /2 I B 2 , j 

PT = 2-^ E E lp±. (34) 

The sum over \x in Eq. (J53)l yields easily 

E + I) 2 = (2I B + 1)(1 + \i B {Ib + 1)), (35) 

and after substitution the value w(Ib) of Eq. (J2H), we arrive at sought for 
result 

)-N B n b/2 



n-i\ B "Bl* / at _i_ 1 \ 4 



ii 



The remaining sum over Ib in Eq. (J36|) depends on whether N B is an even 
or odd number. If N B is an even number, then I™" = and straightforward 
summing over Is in Pi of Eq. (}3l)j) by exploiting the following sums of the 
binomial coefficients, 

N B + 1 \ 




z B =o xl N B + Ib + K 



yields the polarization 

for an odd iV = N B + 1 spin cluster [TT]. If N B is an odd number, then 
I B m = | and some simple algebra gives the polarization 

Pi = m K/) (39) 

for an even N = N B + 1 spin cluster. 

When N > 1 , the T\ of Eq. §^ behaves as (N+2-2(ttN/2)- 1 / 2 )/(3N) . 
Eq-s. (jUBl), (EH give sought for time-independent contributions Pi to the to- 
tal polarization P\{t) of Eq. (|3T)|) for odd and even numbered spin clusters, 
respectively. 

It remains to find the time-dependent contribution P 1 osc (r) to the total 
polarization Pi(r) of Eq. (|2~%j). Among the four pairs (/,/') in Eq. (|2^j). 
only the pairs (J, I') with I ^ I' contribute to the time-dependent part 
of the function Pl(t) of Eq. (j2*S|) . This occurs for I B > | only, since 
otherwise, i.e. for I B = , the allowed values I = I 1 = ~ are encountered 
already in the time-independent polarization Pi . Thus, among the four 
pairs (1,1') in Eq. (|2U|) only the two pairs, i.e. (I = I B + |,P — Ib — \) 
and (I — I B — |,P = Ib + \) are allowed and provide complex conjugate 
contributions to the real- valued function P 1 osc (r) . It suffices to deal with the 
first pair, (/ = I B + |, I' — I B — |) . The polarization becomes 

N B /2 Ib + \ Ib~\ 

p osc (t) = 2 -(iv B -i) £ w{Ib) £ Wx 



-(/ s+ I)m'=-(/ B -|) 
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x2cos(2r(/ i? + i))( V m A cl B - h f , c[ B+ ^ f. (40) 

v v -d 2 V ^— ' ^,mA;lB,m'—mA ^,mA\lB,m—mAJ K 1 

m A =±\ 



To complete the derivation of the function P 1 osc (r) , we use the expression for 
the factor w(Ib) of Eq. (j2U), the C-G coefficients of Eq. (|2^j). and to sum 
in Eq. (|4U|) over the variables m and m' for the fixed value of Ib , yielding 

2 -N B +3 N B /2 / N 1 x x 

- 5(iv^T) £ ( |JVfl : \„ + J W> + 1) + 5)). (41) 

Finally, by gathering the expressions for the Pi in Eq-s. ([38)1 . (j3T?j) and the 
expression P° sc (r) in Eq. (|4*T|) we'll go over (with the substitution k — Ib — 2 
for even TV and k = Ib for odd iV ) to the total polarization on the first 
spin 

A(r) = m Km) + W £ 4* W cos(r(iV - 2*)) (42) 

for an even iV -cluster, and 

JV-l 

TV 4- 2 2 4_iV 2 
A(r) = + — £ 4k(*0 c ^(^(^ - 2*)) (43) 

for an odd iV -cluster, with coefficient 

^. ( ^±i-*)(^-*)(g 

holding in both cases. The formulas (|42jl . (f4*3*|) are the central result of the 
paper, and they can accurately describe a variety of systems in the next 
section. 



4 Discussion 



4.1 Non-Ergodic Spin Dynamics 

As Eq. fj4*3"|) states, for large odd iV -clusters, the time average polariza- 
tion (Pi(r)) of the spin 1 tends to | , while the time average polarization 
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^ of any other spin tends to , i.e. polarization of the spin 1 does not 
spread uniformly over the iV-spin cluster. We call this behavior as non- 
ergodic spin dynamics to confront it with the ergodic spin dynamics pro- 
viding 1/N polarization for all spins in the N -spin ensemble. Figure 2(A) 
shows the behavior of the polarization P\(r) for the series of odd N -clusters. 
The principle features of the periodic pulses of the polarization are managed 
by two factors: the time reversibility of the dynamics affects the exact re- 
entrance of the polarization to the prepared value -Pi(O) = 1 after each 
period Air/g, and, the second, gives rise to the temporary interval with the 
time-independent polarization of the spin 1 , see Appendix for details. For 
large N -clusters, the total period 47r/g can be partitioned into the switching 
time t sw = 4tt^§ and the stopping time t st = y (1 - y^) ; recall r = \gt . 
By referring to the Appendix, the polarization P\(r) is peaked at the mo- 
ments t = 0,2Tc/g, 4:71 /g, .... The profile of the function P\(t) around the 
moment r = is 

Pi(r) = I + |(l ~ r 2 N) exp(-r 2 iV/2). (44) 

The same profile of the function P\(r) holds around all the moments r = 
mix , for all integer m. The interval between the successive peaks and their 
width are 

T=— and A T = 47r-^$i, (45) 
9 gVN 

respectively. In other words, for large N -clusters, N 1 , almost all the 
time the polarization of the spin 1 stays at the fixed value Pi = 1/3. The 
oscillating part of the function P\(r) is the odd function of the time with 
respect to the moments r = |, . . . , as it is apparent from the P\(t) of 
Eq. (USD- 

Figure 2(B) shows the profiles of the polarization for even iV-spin clus- 
ters. For large even values of iV , the polarization at spin 1 stays fixed over 
the long time interval t st = ^(1 — — within each period 2n/g. Unlike 
odd iV -clusters, the profiles P\(r) for even iV -clusters are even functions 
of the time with respect to the time moments r = ir, 2tt, .... 

Relying on the experimental time interval T and the width of the pulses 
At , the two relations in Eq. ()45|) together with the expressions for the 
coupling g of Eq.-s (fi"0|) . (fl~5|) give the volume and form-factor (the angle a 
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is assumed to be known) 

V =-(^Y and /0^(cosa) = -^- 7 , (46) 

where c = N/V denote the concentration of the molecules carrying the 
spin- | . 

4.2 Polarization Dynamics in Fluctuating Nano-bubbles 

The Eq.-s (J4*2~j) . (jlHj) can be adopted to accounting for a time dependence of 
the volume of the nano-cavity and, thus, providing a means to explore NMR 
imaging of the cavitation bubbles in water [20], blood pT] . etc., alongside 
with the conventional high-speed photography. Dynamics of the surface of 
a typical bubbling behavior occurs on a millisecond time scale [20 , i.e. on 
the same time scale which is relevant for the nuclear spin dynamics. It is 
legitimate therefore to question here: how does the dynamics of the nano- 
sized volume affect the nuclear spin dynamics ? Our intention in this section 
is to show that the fluctuations of the nano-volume (governed either by an 
external inputs or by inherently thermal noise) drive the polarization to the 
non-ergodic value 1/3 irreversibly so that the periodic in time pulsating of 
the polarization break down as the time proceeds. 

The formulation of Section El is easily extended to the case of time varying 
volume V since the coupling g(V(t)) enters into the Hamiltonian (|16|) as a 
common factor in front of the whole operator part. The functional form of the 
polarization P 1 osc (r) of Eq. (j4*Tj) which has been derived for time independent 
coupling g is generalized to the case of the function g(t) provided that the 
time t = gt/2 in Eq. (}4~Tj) is now replaced with a new time 

H^/o'^*' < 47) 

We will be interested in the transformation (|4*Tj) 

g(t) = (g) + 6g{t), (48) 

where the Sg(t) stands for the Gaussian random noise characterized by the 
first two moments 

(5g(t)) = 0, {Sg^Sgfo)) = ((5g) 2 ) 7 (|ti - t 2 \), (49) 
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where ((5g) 2 ) is the variance and the j(t) denotes the correlation function, 
for example, j(t) = exp(—t/t c ) , with t c being the correlation time. By the 
comment before Eq. (JUJ), we replace the factor cos \ 2t(Ib + |)J in Eq. (J4*T]) 

with expression cos (i# + \) ({g)t + Jo Sg(t')dt'^ . The Gaussian averaging 
of this factor over the random function 8g(t) is carried out by the formula, 
see e.g. |T5J, 

'exp (i{I B + I) J*8g(t')dty = exp(-(/ B + \) 2 {{5gf) T 2 ), (50) 



with 



T 2 = (\t-t')Mt')dt'. (51) 
Jo 

We confine our attention, first, with the polarization for even N , with N ^> 
1 , and then close the section with the final result for odd N , N ^> 1 . Let 
us write the polarization of Eq. (|41j) with the averaging Eq. (J5U)) as 

■z 

<&(*) = exp(-(/ B + i) 2 ((^) 2 )T 2 )cos((^(/ i? + i)). (52) 




The exponent in Eq. ()52|) tell us that the successive peaks of time dependent 
part of the function P\ (t) , reduce to zero at t —>■ oo , thus, only the time 
independent part, i.e. Pi = 1/3 of Eq. (f3T)J) of the function Pi(t) survives at 
t — ► oo , after the Gaussian averaging over the function 5g(t) . The integral 
over t' in the constant T 2 of Eq. (JoTj) can be evaluated in the two asymptotic 
cases, for the large and the small correlation time [To] . 

2 _/t 2 /2, if t c 2 ((^) 2 )»l 
~W, if t c 2 ((^) 2 )«l ' [b6) 

In order to find the function P\(t) of Eq. (|52j) at asymptotics N ^> 1 , we 
can replace the sum in Eq. (J52)) with the Gaussian averaging, just as done 
in Eq. (A. 3) of the Appendix, yielding 

1 16 ^ 2 „ , , / , 1 



Pi(t) = tt + -7= E ™s((g)tn) n 2 - - )e~™\ (54) 
3 3^/2^/2 t^i V 47 
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where 



| + m 2 )T\ (55) 



Figure 3 shows the polarization dynamics of a single spin within N = 134 
spin aggregate for ((6g) 2 )/ (g) 2 = 1CT 4 . Based upon the Appendix, the sum 
over n in Eq. (jSljl is simplified with the Poisson's summation formula. On 
defining the partial sums entering into Eq. (|54p by 

N/2 

S x {t) = Y, c os((g)tn)e- an2 = 

71=1 

= -5 + 5^1? M-Vf )'/<■). < 56 > 



and 



N/2 , 

S 2 (t) = £ cos((^)to)n 2 e- an = -dSi(t)/da = 

n=l 

s£r Jjl - + ^) 2 ) exp(-^( 9 4- ^) 2 /a), (57) 



we get 



PlW ^ + S^( S2(() -3 Sl(() )- (58) 

To find the envelope of the successive peaks of the function P\(t) of Eq. 
(fESjl. put into Eq.-s (jo^j) - (J58J) the time t = 27im/(g) , with m running over 
integer numbers. This gives the polarization at the discrete moments m , 

N / 2 1 1 oo 

ft(m) = £ ^ = 4 + l£ E e " 7r2fc2/Q ' (59) 

71=1 fc = — OO 



&(m) = £ « 2 ^" 2 = T^J E (1 " (60) 

n=l k= — oo 

The functions Si(m) and 5*2 (to) inherit their dependence on the "time" m 
through the constant a of Eq.-s fl53j) . (jUj) 

f| + 2vr 2 m 2 ((^) 2 )/^) 2 , if t c 2 ((^) 2 )»l 
\f + 2 7 rm* c <(^) 2 )/<< ? >, if t c 2 ((W)«l ' 

17 



where we substitute t = 2iTm/{g) into Eq.-s ([55)1. (j5HJ). For N ^ 1 and 
m ^> 1 , we drop the summand 2/N in Eq. (J61j) assuming that a ^> 1 , 

r27r 2 m 2 ((W)/^) 2 , if t c 2 ((W)»l (m] 
a ~\2nmt c ((5g) 2 }/(g}, if t 2 ((<^) 2 > « 1 " ^ 

For a ^> 1 , we find the sums over the in Eq. (jH^j) and Eq. (jHOj) by making 
use, again, the Poisson's summation formula, Eq. (A. 5) of the Appendix, 
that accelerates the convergence of the sums for a ^> 1 . On reading Eq. 
(A. 5) of the Appendix backward from the r.-h. side to the l.-h. side, we get 



1(a) = J- e- n2k2/a = e~ ae2 = l + 2e- a + (D(e~ 4a ). (63) 

' 12 k=~oo l=-oc 

Thus, Eq.-s (J5fi|) . (JoTj) become 

Si(m) = -\ + \l{p) = e~ a , (64) 

so that Eq. ()58|) gives the polarization of the first spin 

with the a from Eq. (|62j) . Eq. (JfiHj) implies the total polarization of the first 
spin 

We conclude this section with the result for the total polarization for odd 
total number, N , of spins. Due to alternating peaks of the polarization 
P\(t) of Eq. (J4*3~j) . see also Figure 2A, we get 

!l — ^T7^7= e~ a i for large odd m 
3 + jv3/2^: e > tor large even m 

Eq. (J58j) with the a from Eq. shows that the the polarization peaks, 
Pi(m) , of spin-carrying gas has the Gaussian and the exponential time de- 
pendence for the large and the small correlation times of the fluctuations of 
the nano-bubbles, respectively. 
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4.3 NMR Line Shape 

To calculate the NMR line shape exactly, we use the same effective Hamil- 
tonian of Eq. (JTfij) as described in Section EJ The NMR line shape is the 
Fourier transform of the free induction decay (FID), F(t) , of the iV-spin 
ensemble jT^]. It is the NMR line shape on the protons in hydrogenated thin 
silicon films that got the first experimental evidence for the validity of the 
effective Hamiltonian of Eq. ()16|) in the nano-cavities |4J. We are interested 
in the FID signal 



F(t) = tr {e tHt I+e- im I.}/tr {/+/_} , 



(69) 



With I a = En=l I 



y , ex x,y,z 



The reason for an exact 



solution for the FID of Eq. (JfiTlj) is that the total Hamiltonian of Eq. (j!6|) can 
be expressed in terms of the 3 collective spin operators I a , instead of 3N 
spin operators in H of Eq. (fi7)|) . Due to commutation relations [I 2 , /J =0, 
rewrite Eq. (JBTJj) as 



F{t) = tr {e iGtI * 



F}/tr {/+/_} 



(70) 



with G = 3g/2 for dipolar interactions in the effective Hamiltonian of Eq. 
()16|) with ( = 2. The Heisenberg equations of motion for operator /+(£) = 
" is solved exactly to I+(t) = exp(iGt(2I z - l))/+(0) , I+(0) = 



AGtii 



iGtii 



I + . The averaging in Eq. 
basis of the states 
N l = N and I z = (iV T 



N h N x 



is performed in N\/(Nf\N^\) -fold degenerate 
with iV-f , (Ni) spins up (down), so that jVj + 



Ni)/2. The averaging gives the FID 
3 



F(t) 



cos 



:9t 



N-l 



(71) 



The effect of dephasing of the proton spins within the nano-cavity due to the 
interactions with the protons at the surface of the nano-cavity is introduced 
phenomenologically as 



F(t) 



cos 



3 

2' jt 



N-l 



-t/T 2 



where, the time T 2 relevant for the experiments [4 is T 2 
moments of the line shape are 



duu n %(u) 



d n F(t) 
d(it) n 



(72) 

1-^3 ms. The 
(73) 



t=Q 
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where enters through the Fourier transformation of the FID, 

/oo 
du^{u)e lwt . (74) 
-oo 

The second and the fourth order moments are inferred from Eq. (J71|) . as 

M 2 = (N- 1) (3g/2) 2 , M 4 = (N - 1)(3JV - 5) (3^/2) 4 . (75) 

The M2 derived in j3] by the Van-Fleck formula coincides with the M2 of Eq. 
(J75|) . as it should be. The line shape in the nano-cavity volume appears to be 
the volume dependent allowing for determining the volume of the nano-pores 
in hydrogenated silicon films |3]. 

5 Conclusion 

We have presented the exact time-dependent description of spin- \ dynamics 
with infinite range spin interactions and the initial polarization prepared on 
a single spin 1 , i.e. -Pi(O) = 1 . Spin dynamics for odd and even numbered 
clusters demonstrates the periodic pulses of the polarization P\{j) on the 
spin 1 . For large odd N -clusters, the polarization on the spin 1 has pulses 
over the time interval t sw = 47r^= , from -Pi(O) = 1 to the time-independent 

polarization which lasts, thus, t st = ^(1 — ^7^) within any period 47r/g. 
For large even N -clusters, the switching time t sw = 27r^== , and the period 
equals 27r/g . The stationary polarization on the spin 1 is non-ergodic, as 
its value tends to 1/3 (instead of tending to the ergodic value 1/N) as N 
tends to infinity. The profiles of the polarizations within the series of odd 
(even) large clusters are remarkably similar. 

The specific polarization profile in the clusters with infinite range spin 
interactions is in a sharp contrast to the profiles of the polarization in ID 
clusters with nearest neighbor XY Hamiltonian [3]. Two differences can be 
drawn from the presented results: 

(1) The overall behavior of the polarization Pi(t) in the system with 
infinite range interaction is strictly reversible, periodic with the period 47r/g 
for any iV , whereas, on a large ID chains, N 1 , with XY Hamiltonian 
the polarization P\{t) on the spin 1 moves in irregular fashion. 

(2) For large iV-spin clusters, iV ^> 1 , the polarization Pi(t) of the spin 
1 exhibits the plateau region at non-ergodic value Pi = 1/3 and the pulses 
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of the polarization Pi(t) have a short time span of about ATTO(l)/(gy/N) . 
This is in contrast to the behavior of the polarization Pi(t) in ID spin 
chains with XY Hamiltonian where polarization on the spin 1 depends on 
the time in irregular fashion with t st = . 

Finally, the paper demonstrates the sensitivity of the polarization dy- 
namics ( reversibility and ergodicity in the many-spin systems) to the radius 
of interaction. Incorporation into the theory of the real dipolar interactions 
is the most challenging task of dynamical theory and the accurate answer is 
not settled yet, although the general picture of the spin dynamics is known 
to be diffusional |22j. 
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Appendix. Derivation of Eq. (I44|) 

We want to prove that the function P° sc (t) in Eq. (|4"T|) for N ^> 1 has 
the form of pulses of width 



gVN 

at equidistant time moments r = 0, 2tc, Atc, ... so that the profile of the 
function P° sc (r) at the time moment r = is 

P° sc (r) = -(l - r 2 iV) exp(-r 2 iV/2), for iV > 1 (Al) 

To prove Eq. (A.l), we introduce a new variable n = Ib j t\ i n the expression 
P° sc (r) of Eq. (gU), so that the function P° sc (t) takes the form ( recall that 
the total number of spins equals A" = Nb + 1 ) 

^W^| 2 -(/ + J(" 2 -i)co S , 2 ™,. (A2, 
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Next, we use asymptotic formula for the binomial coefficient, 



■«( » \- 1 exp (--£- (A3) 



jA + n, 



^JttN/2 




Eq. (A. 3) allows one to look at the summation in Eq. (A2) as an averaging 
over the Gaussian distribution function. To simplify the calculations of Eq. 
(A2) further, we apply the Poisson's identity 



£=— oo ' ^ k=— oo 



]T cos(2ne£)e- ae2 = J- £ e - 2 (^)> ( A4 ) 



In many circumstances, including our ones, the resulting sum over the k in 
the r.-h. side of Eq. (A. 4) converges much faster than the original sum over 
the t in the l.-h. side of Eq. (A.4). To apply Eq.-s (A.4) to Eq. (A. 2), we 
expand the sum in Eq. (A. 2) up to n = oo since the terms of the sum in Eq. 
(A. 2) practically vanish for n > N/2 and N ^> 1 . Thus, by the Poisson's 
identity of Eq. (A.4), we introduce the sum ( a partial contribution to the 
sum of Eq. (A. 2) ) 



Sl (r) = £ cos (2rn)e-^ = - 1 - + 1 - ] S ^ £ «p(-V(* + 

n=l ' k=—oo 

(A5) 

In order to show that the function Si(t) has the form of the Gaussian peaks 
at equidistant moments r = 0, ±7r, ±2ir, ... it suffices to analyze the function 
Si(t) around the point r = . In this case, the leading contribution to the 
sum in Eq. (A. 5) is provided by the term k = . Denote that if we analyze 
the peak around r = inn , where m is integer, then the leading contribution 
to the sum S\(t) in Eq. (A. 5) comes from the k = —m. Thus, in considering 
N ^> 1 , we can drop all the terms in Eq. (A.5), except the leading term 
k = , yielding 

S 1 (r) = - 1 -+ 1 -J^e-^. (A6) 



Analogously, we determine the partial sum 

(JV-l)/2 



S 2 (t) = E" cos (2rn) n 2 e~ n 

71=1 

_ a 



5 1 (r) = ^(l-r 2 A)e^/ 2 . 



(A7) 



£>(2/7V) wi V' J sV2 
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At iV 3> 1 , the function Si(t) of Eq. (A. 6) has a negligible contribution to 
the function 

Pr(r) = 1 -^ T = (s 2 (t) - ]s 1 (r)) 

in comparison with the contribution of the function S 2 (r) of Eq. (A. 7), 
yielding the sought for result of Eq. (A.l). 

In general, the function P 1 osc (r) for an arbitrary r has the pulses at the 
moments r = kn with integer k , 

PT(r)= £ ^(l-vr 2 (A; + ^) 2 iV)exp(- 7 r 2 (A; + ^) 2 f). (A8) 
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Captions to figures. 

Fig. 1 . Schematic representation of nano-pore with the molecules carrying 
the nuclear spins and are in a rapid thermal motion. 

Fig. 2 . On the panel A, the polarization Pi{r) of Eq. ()43|) of the first 
spin is varied with non-dimensional time r = gt/2 for the series of an odd 
total number, N , of spins. The panel B displays the polarization P\(t) of 
Eq. (J42j) for the series of even N . 

Fig. 3 . Polarization dynamics of a single spin, Pi(t) , within the N = 134 
spin aggregate when the volume of the nano-cavity fluctuates providing with 
relative variance for the g coupling equal to ((Sg) 2 ) / (g) 2 = 10~ 4 , see Eq.-s 
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spin-i/2 probe 
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